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Abstract 

We report the experimental observation of two-mode squeezing in the oscillation quadratures of 
a thermal micro-oscillator. This effect is obtained by parametric modulation of the optical spring in 
a cavity opto-mechanical system. In addition to stationary variance measurements, we describe the 
dynamic behavior in the regime of pulsed parametric excitation, showing enhanced squeezing effect 
surpassing the stationary 3dB limit. While the present experiment is in the classical regime, our 
technique can be exploited to produce entangled, macroscopic quantum opto-mechanical modes. 

PACS numbers: 42.50.Wk,45.80.-br,05.40.-a,07.10.Cm 
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Cavity opto-mechanics [1] is achieving several major experimental breakthroughs, includ¬ 
ing the optical cooling of opto-mechanical oscillators down to a thermal occupation number 
around or below unity [2-6]. A great deal of attention is now devoted to possible realizations 
of strongly non-classical states of macroscopic variables, that would open the way to new 
quantum information tools, as well as to crucial tests on the classical-to-quantum transition. 
In this framework, several techniques to obtain squeezing of the motion quadratures have 
been studied [7-12] and demonstrated on thermal oscillators [13-19], including quantum 
non-demolition measurements with modulated [7, 15] and pulsed [10, 17] readout, and para¬ 
metric excitation [11, 13, 14, 16, 18, 19]. Very recent reports describe the quantum squeezing 
of the motion of a cooled nano-oscillator exploiting the mechanical interaction with a mi¬ 
crowave held [20, 21]. A different class of macroscopic, strongly non-classical systems can 
be obtained by involving more than one optical and/or mechanical mode. In particular, two 
mechanical modes have been exploited in opto-mechanical hybridization [22] and are the 
basic ingredients of two-mode mechanical squeezing, that is the subject of this work. 

A two-mode mechanical squeezed state is characterized by correlated huctuations be¬ 
tween one quadrature of the hrst oscillating mode and one of the second mode (as well as 
between the corresponding conjugate quadratures). In a quantum system, such correlations 
translate into entanglement between the two oscillators, i.e., the mechanical correspondent 
of optical twin beams generated in optical parametric amplihers [2]. In a suitable combina¬ 
tion of quadratures of the two oscillators, the huctuations are reduced below their standard 
quantum level (or, in a classical system, below the thermal noise level). Such a reduction 
can be obtained by implementing a modulation of the spring constant at the sum of the res¬ 
onance frequencies of the two oscillators. Experimental demonstrations of fully mechanical 
two-mode squeezing of thermal huctuations have been recently reported for nano-oscillators 
where the parametric modulation is obtained by electro-mechanic driving, either incorporat¬ 
ing a piezoelectric transducer directly into the mechanical elements of a device designed on 
purpose[24], or taking advantage from an accidental resonance with a mode of the substrate 
of a membrane oscillator [25]. Here we study two-mode thermal squeezing in a cavity opto¬ 
mechanical system, exploiting parametric modulation of the optical spring. Our method, 
just based on the fundamental opto-mechanical interaction, does not require peculiar oscil¬ 
lator characteristics and can be extended to any cavity opto-mechanical system, including 
microwave driven nano-mechanical oscillators, and even in a moderately resolved sidebands 


2 


regime. In addition to stationary variance and spectrally resolved measurements, we describe 
analyses of dynamic behavior in non-stationary conditions, showing enhanced squeezing ef¬ 
fect surpassing the stationary limit determined by the onset of parametric instability. 

Our silicon oscillator is made on the 70 pm thick device layer of a SOI wafer [26] and is 
composed of a central mass kept by structured beams [27], balanced by four counterweights 
on the beams joints [28]. On the surface of the central disk, a multilayer Si 02 /Ta 205 dielec¬ 
tric coating forms a high reflectivity mirror. For the main oscillator mode, labeled ® in Fig. 
la (resonance frequency 172 kHz, effective mass 250pg), the motion of the mirror is balanced 
by the motion of the counterweights. As a result, the total recoil on the inner frame is null 
and the joints become nodal points. In a second oscillation mode, the counterweights move 
in the same direction of the mirror (mode ® in Fig. la), giving slightly different frequency 
(225 kHz) and mass (lOOpg). Thanks to the overall design strategy, both modes repro- 
ducibly display a mechanical quality factor of about Q ~ 5 x 10^ at room temperature, only 
limited by the thermoelastic dissipation in the substrate, and the balanced mode exhibits a 
Q exceeding 10® at cryogenic temperature [29]. The oscillator reflecting surface is the end 
mirror of a Lcav = 1-27 mm long, high Finesse {T = 18000, half-linewidth k/2'k = 3.2 MHz), 
over-coupled (input mirror transmission T = 300 ppm) Fabry-Perot cavity. 

In our setup (Fig. lb), two laser beams derived from the same Nd:YAG source are 
superimposed with orthogonal polarizations and coupled to the cavity. The hrst one is used 
to obtain a signal proportional to the detuning, with an rf sidebands technique. This signal 
is used to lock the laser to the cavity resonance, with a bandwidth of about 10 kHz, and 
to measure the oscillator displacement x{t). Furthermore, it is sent to two double-phase 
digital lock-in amplihers whose outputs are acquired for the reconstruction of the motion of 
the oscillating modes. The second, more powerful beam, with a controllable frequency shift 
with respect to the previous one, sets and controls the optical spring. 

The intracavity radiation pressure force, depending on the detuning and therefore on 
the oscillator position, originates a so-called optical spring [1] that modihes the mechanical 
resonance frequency. Moreover, due to the hnite cavity held buildup time, it changes the 
effective damping of the oscillator motion, resulting in modal optical cooling (cold damping) 
when the input radiation frequency 0JL/27r is red shifted with respect to a cavity resonance 
u}c/2n. In the bad cavity limit and for small detuning (i.e., when both the opto-mechanical 
frequencies and A = ool — ooc are well below k), the optical spring constant can be approx- 
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FIG. 1: (Color online) (a) Probability distributions in the (Yi,X 2 ) quadratures plane. Left: thermal 
distributions, without parametric excitation. Right: with parametric excitation, the fluctuations 
in the two quadratures are strongly correlated. Vi and X 2 are normalized quadratures of two 
different oscillating modes, whose shape is shown by Finite Elements Model simulations on the 
left of the panels, (b) Experimental setup. The hrst oscillation mode is parametrically locked at 
fo using the optical spring. The frequency fsum-fo) corresponding to the resonance of the second 
mode, is selected after a mixer by the low-pass hlter (LPF) and used as reference for the second 
lock-in ampliher. 
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imated as -ft'opt — where P is the intracavity power. In the same conditions, the 

effective linewidth becomes 7eff = 7m + 7opt with 7opt — — is the mass, uj^/2tt the 

resonance frequency and 7 the linewidth of the free oscillator) [1, 30]. We exploit the optical 
spring to stabilize the frequency of the first oscillation mode, by slightly exciting it with an 
amplitude modulation on the strong laser beam, comparing its motion with the reference 
oscillator on the first lock-in, and phase-locking its oscillation to the reference by acting on 
the optical spring through the pump laser detuning. Such parametric locking [18] is also 
indirectly efficient on the second mode, since it stabilizes the whole optomechanical system. 

We now consider two oscillation modes, both changing the cavity length, with effective 
displacement Xi{t) {i = 1,2), momentum Pi{t), resonance frequency cuj and mass m*. It is 
convenient to observe the oscillators in their respective rotating frames. We thus analyze 
the behavior of the slowly varying quadratures, defined as Xj = Xi cos cujf — sin ooit and 

Yi = Xisinuiit -|- ^^^coscujf. We remark that the parametric locking above described acts 
on one quadrature of the first oscillator (namely, Xi), while the orthogonal one remains 
unaffected. Therefore, in the following we just compare the model with the experimental 
data concerning the Yi quadrature, even if for the theory (neglecting active parametric 
locking) the two quadratures are formally equivalent. 

In our experiment, the laser detuning is modulated at oji + U 2 such that A —)■ A[1 -|- 
(3 cos(a;i + uj 2 )t] (with (3 -C 1). As a consequence, the optical spring constant is also modu¬ 
lated, and the resulting force is —(xi -|- X 2 ) K/s cos(ci;i -|- u) 2 )t, with = Xopt/d- 

Writing the modulated force in terms of the quadratures, we derive that its quasi-resonant 
terms are associated to the two-mode squeezing Hamiltonian Hj = ^ (X 1 X 2 — hiY 2 ), that 
couples the quadrature Xi with Y 2 and hi with X 2 . As a consequence, the two-mode system 
can be decomposed into two equivalent independent subsystems. The Langevin equations 
governing their motion can be written as 


X + AX = £ 


( 1 ) 


where X = (X, F) can be either (Xi,F 2 ) or (yi,X 2 ), ^ with 

{fi{t)fj{t')) = 5ij5{t — t'), and Spi is the force spectral density acting 

on the f-th oscillation mode. The system matrix is 


A = 




( 2 ) 
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with gi = Kp/Auirrii. 

Stationary, null-average solutions are stable if ^ = Kj^/Kth < 1 where the threshold Kth, 
calculated by nulling the determinant of A, is Kth = 2^'ji'j2'mim2UJiU]2- In this situation, it 
is interesting to evaluate the peak spectral densities of the quadratures (at cn = 0). Indeed, 
they are proportional to the variance measured after integration over periods much longer 
than 1 / 7 , that is a frequent experimental procedure. 

In the simplihed case of oscillators with identical susceptibilities (with mass m, resonance 
frequency Uq, damping rate 7 , force spectral density Sp), the eigenvalues of A are simplified 
to A± = |(1 ±^) and the respective eingenvectors are X± = (A ± Y)/\/2. Below threshold, 
the spectra of the fluctuations around the stationary solutions are 


c ( \ 2Aj- 2 


(3) 


where the variances are = (cr°)^ 7 / 2 A± with (a^)^ = Sp/{2'ym?u:Q). On the combined 
quadrature A+, the peak spectral density S'x+(0) is reduced by a factor of (1 -|- g^. The 
total variance is accordingly reduced by a factor of (1 + g), with a maximum reduction 
of respectively - 6 dB and -3dB when approaching the threshold. The quadrature X_ has 
maximal fluctuations, with peak spectral density and variance increased respectively by 
1 /(1 —and 1/(1 —^) with respect to their values in the absence of parametric modulation. 

In the realistic case of different parameters between the two oscillation modes, to main¬ 
tain meaningful combined quadratures it is useful to first normalize each quadrature to its 
standard value, calculated without parametric excitation. Therefore, we calculate the nor¬ 
malized quadratures Aj oc Aj and F* oc Yi in order to have either S'- (0) = S'- (0) = 1 

-Aj li 

or (a - = (cr~ = 1 , according to the chosen indicator (we use the ” 0 ” superscript for 
the quantities measured without parametric excitation, i.e., at = 0). Eqs. (8-10) remain 
formally the same, with a different physical meaning of the parameters gi and S'^. [31]. The 
general combined quadrature is now defined as A^ = A cos 6 *-|-F sin^, with spectral density 
Se and total variance aj. Without parametric excitation, we have symmetric distributions 
with either 5^(0) = 1, or {aj)"^ = 1, for any value of 6, as shown in the left panel of Fig. 
la. With parametric excitation, the two quadratures corresponding to eigenvectors of A 
still maintain a Lorentzian spectrum, with width A±. However, differently from the case 
of identical oscillators, they do not correspond to extrema in the quadrature fluctuations. 
Either S'e(O) or aj should be minimized with respect to 9 to derive the combined quadrature 
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exhibiting the largest squeezing, in correspondence of 6'min- The probability distribution 
with parametric driving is shown in the right panel of Fig. la. The strong anti-correlated 
fluctuations in the two quadratures are clearly revealed by the elongated shape of the graph. 
In a narrow interval of 6, close to vr/d, the fluctuations of Xg are below the thermal ones. 

It is important to remark that the described picture assumes an exact detection phase. In 
the realistic case of (small) fluctuations A6 during the measurement time, the measured min¬ 
imal spectral density is + ((A6')2)S'max where S^m = Sg^,^{0) and = Se^,^+n/ 2 { 0 ) 
(a similar expression holds for the variance). This feature is crucial to describe the experi¬ 
mental results. The expected minimal spectral density now reaches a minimal value before 
the threshold gain, followed by a deterioration of the squeezing. 



0.0 0.2 0.4 0.6 0.8 1.0 

Parametric gain g 


FIG. 2: (Color online) Experimental measurements of the spectral densities 5'min (full symbols) 
and Fmax (empty symbols) as a function of the parametric gain, together with the corresponding 
theoretical predictions (solid lines). The signal from the four quadratures are preliminarily nor¬ 
malized in order to have unity spectral density in the absence of parametric excitation. Different 
colours (grey tones) corresponds to two data sets, taken with different levels of optical cooling. A 
dash-dotted line shows the prediction for a perfectly stable and optimized detection phase. 

In Fig. 2 we show our experimental measurement of S'min and S'max as a function of the 
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parametric gain, together with the corresponding theoretical predictions. The experimental 
values of S'min and 6 'min are obtained by globally minimizing the peak spectral density of the 
combined quadrature {Xi cos 0i + Yi sin 0i) cos 6 + (X 2 sin 02 + ^2 cos ^ 2 ) sin 9 with respect 
to 01 , 02 and 6, to account for possible errors in the hxed lock-in reference phases. 

As expected and already observed [24, 25], the squeezing does not surpass - 6 dB before 
the onset of parametric instability. Moreover, S'min reaches a minimum before the threshold 
gain, in contrast with the expected behaviour for a perfect and stable detection phase, dis¬ 
played in the picture with a grey dash-dotted line. This feature, already observed in recent 
experiments [24], is well reproduced by the complete model with the parameters of our differ¬ 
ent oscillating modes. The theoretical values are calculated using independently measured 
oscillators parameters, and just htting to the experimental data the threshold parametric 
gain, whose independent estimation is not enough accurate, and the phase variance ((A 6 *)^) 
that results to be around 10 “^ rad^. 

Above threshold, the system admits no stable solutions and therefore no stationary spec¬ 
tra. However, one can still calculate the behavior of the variance in measurements of a 
quadrature, performed at the time t after the switch-on of the parametric modulation [31]. 

To start describing the physics involved, we present again the case of identical oscillators. 
As before, the extrema of the variance are found for 6 = ± 7 r /4 at any time t, i.e., in 
correspondence of the eigenvectors of A, with and = g^. The respective 

temporal evolutions are g±(t) = (g°)^[exp(—2A±t) -|- 2^(1 “ exp(—2A±t))]. For ^ > 1, 
Umin approaches 1/(1 -|- ^) (i.e., a squeezing exceeding -3dB) with a time constant Tshort = 
1/2A+ = 1 / 7(1 +g), while exponentially diverges with the longer time constant pong = 
1/2|A_| = /'^{g — 1), exhibiting the critical slowing down at threshold. If the system 

dynamic range is large enough, it can be left evolving after the switch-on of the parametric 
excitation for few Tshort, reaching therefore a stronger two-mode squeezing with respect to 
the stationary situation. 

For different oscillators, the variances a± in the combined quadratures corresponding to 
eigenvectors of A still evolve exponentially with time constants 1/2A±, while for any other 
quadrature we expect the combination of exponentials with time constants 1/2A+, 1/2A_, 
and 1/(A+ -|- A_). As in the stationary case, minimal fluctuations are no more found for the 
eigenvector corresponding to A+ (except for asymptotically long t, as we will see). Moreover, 
here the optimal phase 9 depends on t, as well as on g. 



FIG. 3: (Color online) Time evolution of the experimental minimal variance after the onset 
of the parametric modulation (symbols), together with the prediction of the model (solid lines), 
for different modulation amplitudes. The normalized parametric gains g extracted from the fits 
are g = 0.92 (a), g = 1.06 (b), g = 1.30 (c). The second fit parameter is the phase noise, that 
results to be A0rms = 0.041 rad (a), A^rms = 0.045 rad (b), A^rms = 0.070 rad (c). The variance 
is measured over 5000 consecutive bursts, repeated at the rate of 50 Hz. The duty cycle is 50% for 
panels (a) and (b), and 15% for panel (c). Light (green) dashed lines: predicted evolution in the 
absence of phase noise {A6 = 0). Black dashed lines: predicted evolution of 

We have applied to our system periodic bursts of parametric modulation. The bursts are 
long enough to observe the system evolution, and are repeated at the rate of 50 Hz. The 
four output signals of the two lock-in amplihers (integrated with a time constant of 40 /rs) 
are sampled at 1/At = 200 kS/s, then sectioned synchronously with the periodic bursts, so 
that each section starts at the beginning of a burst. We note that the sampling rate is much 
larger than A±. For each n-th data sample, we calculate variances and cross-correlations 
between the four signals, at the time nAt after the onset of the parametric modulation, 
over 5000 consecutive bursts. The minimal variance is obtained by global minimization 
with respect to the three phases 0i, 02 and 6. In Fig. 3 we report the time evolution of 
crTn, together with the prediction of the model, for three levels of gain g below, near and 
above threshold. As expected, in the last case the variance falls below the stationary limit 
of ~ 0.5. The effect of the phase fluctuations A6 is even more crucial than in stationary 
conditions. The measurable minimal variance presents a minimum for an optimal value of 
t, after which it starts increasing due to the mixing with the exploding quadrature. We also 
remark that, in the absence of phase fluctuations, there is a qualitative difference between 
the two regimes ^ < 1 and g > 1. In the former case, the minimal variance is always lower 
than cr^, as in the stationary case (see Fig. 3). Above threshold, instead, it asymptotically 
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converges to and the combined quadrature with minimal fluctuations tends to coincide 
with the stable eigenvector. 

In conclusion, we have prepared an optomechanical system in a two-mode squeezed ther¬ 
mal state. Some features recently observed in purely mechanical and electro-mechanical 
systems, such as the deterioration of the squeezing before the parametric threshold, are ac¬ 
curately observed and well reproduced by a model considering different mode parameters. 
Moreover, we analyze the dynamical evolution after the onset of the parametric modulation 
and show that the -3dB limit in the squeezing of the two-mode variance, rigorous in sta¬ 
tionary conditions for identical oscillators, is overcome for a bit of time. Even in this case, 
our model well reproduces the experimental data, and underlines the crucial role of phase 
fluctuations. Indeed, both the achievable minimal uncertainty and its temporal stability are 
actually determined by such phase noise. 

Our system operates with low optical and mechanical losses and implements optical 
cooling and frequency control, with a setup similar to those employed to prepare the oscillator 
close to its ground state. In such conditions, our technique could produce entangled states 
of the mechanical oscillation modes [32, 33], an intriguing result in the held of quantum 
macroscopic systems. The model here presented neglects the held fluctuations, privileging 
a clearer focus on the main features that are demonstrated experimentally. However, a 
complete quantum analysis shows that peculiar quantum effects can be created involving 
the electromagnetic held, efficiently coupled to the cavity and detected. Thanks to its 
flexibility, our method opens the way to more elaborated schemes including more than two 
oscillating modes, overcoming the experimentally challenging requirement of equal oscillation 
frequencies, and leading to experimental studies of complex multipartite entanglement in 
macroscopic systems. The experimental techniques here described can also be applied to 
different devices embedded in the same optical cavity, or even in separated cavities sharing 
the same optical field [34], thus producing well separated, entangled macroscopic oscillators, 
particularly interesting for investigating quantum decoherence. 
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and by INFN (HUMOR project). A.B. acknowledges support from the MIUR under the 
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SUPPLEMENTAL MATERIAL - DYNAMICAL TWO-MODE SQUEEZING OF 
THERMAL FLUCTUATIONS IN A CAVITY OPTO-MECHANICAL SYSTEM 

The system of two different oscillators 

The quantum Langevin equations for the position x{t) and the momentum p{t) of a 
mechanical oscillator with mass m, damping rate 7 and resonance frequency Uq are 


X — p/m = 0 
p + mcug X + my x = 


( 4 ) 

( 5 ) 


with {f{t)f(t')) = 6{t — t') (valid in the limit of high mechanical quality factors Q [1]). The 



is the thermal mean occupation number. In the limit of high temperature {k-^T hwo), 


valid for the present experiment, we obtain the classical thermal noise spectral density Sp 


2kBTmooQ/Q, while in the opposite case of low occupation number (n 1 ) the quantum 
noise is Sp — 2hmul/Q. If the oscillator is embedded in an optical cavity, both its resonance 
frequency and its damping rate are modihed by the opto-mechanical interaction, therefore 


12 



ojq and 7 must be intended as effective parameters of the opto-mechanical system. The most 
meaningful effect is on the damping rate (that is usually strongly modihed), therefore we 
keep distinguished the oscillator effective resonance linewidth 7 from its natural linewidth 
ojq/Q. In the classical limit, valid in our experiment, the additional noise source due to the 
quantum fluctuations of the radiation can be neglected, and the oscillator can be considered 
at thermal equilibrium with an effective temperature Tg// = Tuq/'jQ. When the quantum 
regime is approached by optical cooling (Teff approaching hwo/kB), held fluctuations become 
important and a limit effective temperature is achieved. Our simplihed calculations remain 
meaningful by replacing ^ ^ 7 and T T^ff. 

It is useful to describe the motion of the oscillator using the slowly-varying quadratures 

X and Y, dehned as X = x cos unt -^sinwnt and Y = xsinwnt -I — —cosuot. Their 

evolution equations (neglecting the terms rotating at ~ 2uQt) are 


2 V 2 mujQ 

(6) 

Y= ^Y+\l^ ^ 

2 V 2 muo 

(7) 


with = 

We now consider the system of two different oscillators (singled out by the subscripts 
”1” and ”2”), coupled by a force of the form —(xi -t- X 2 ) Kjs cos(a;i -|- u) 2 )t. Among the 
different force terms that are obtained by writing xi and X 2 in terms of the respective 
quadratures and expanding the periodic functions, the relevant (quasi-resonant) ones are 
(F 2 sina;it — X 2 Cosa;it) for the hrst oscillator and (Yisina; 2 t — Xicosa; 2 t) for the 
second one. The evolution equations for the quadratures, neglecting again the force terms 
oscillating at twice the resonance frequency, become 

(9. + y) 

(a + y) n + 

mjuj {d, + y) X2 + 

" 12^2 (d, + y) yi + ^Xi = 

where the stochastic force sources have correlation functions = 6ij6abS{t — t'). 
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The coupling terms are associated to the Hamiltonian Hj = ^ {X 1 X 2 — YiY 2 ). Introduc¬ 
ing the ladder operators in the rotating frame as Oj = we can also write 

the two-mode squeezing Hamiltonian in the form Hi = ^0102 + ^ 1 ^ 2 ^ [2] where 

g = Kp/ 2 yJ''yi'j 2 mim 2 UJiUJ 2 is the parametric gain normalized to threshold. Hj couples the 
quadrature Xi with Y 2 and Yi with X 2 . As a consequence, the two-mode system can be de¬ 
composed into two equivalent independent sub-systems. The Langevin equations governing 
the motion of each sub-system can be written as 


X + AX = £ 

where X = {X,Y) can be either (Xi,l 2 ) or (Yi,X 2 ), 




2 mii^i • 

Sf2 1 


2 m2UJ2 • 

with {fi(t)fj(t')) = 6 ij 6 (t — t') and the system matrix is 

71 


A = 


9i 

n T'2 

92 -2 


with gi = The eigenvalues of A are 


A.uJimi 


A+ = 


7i + 72 


± 


7i - 72 


+ 7i72 S' 


and they are simplihed to A± = |(1 ± g) for identical oscillators. 


( 8 ) 

( 9 ) 


( 10 ) 


( 11 ) 


Stationary solutions 


Below threshold, we can calculate stationary spectral densities and variances for any linear 
combination of the quadratures X and Y. We start by calculating the spectral densities at 
vanishing frequency. At this purpose, we neglect the hrst term on the left of Eq. (8). 
We normalize the quadratures in order to have unitary spectral density in the absence of 
parametric excitation (i.e., for gi = 0). We dehne therefore X = X/^/S^ and Y = Y/^/^ 
with = 2S'i7’i/(mia;i7i)^ and Sy = 28 ^ 2 !{''^^ 2 ^ 212 )'^■ Eq. (8) can be written in the simple 


form 


( I ga 
yg/a 1 




( 12 ) 
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where we have dehned a dimensionless asymmetry parameter a = classi- 

^ ^ ir Y SFl'm2CJ272 

cal (i.e., high temperature) limit, the asymmetry parameter can be written as \/ 

while in the strongly quantum regime (with ^ 7 * and n 1), a tends to 1 even for 

different oscillators. The solution of Eq. 12 is 


and has power spectral densities 


Y 


fl-ga /2 

l-§2 

/2-g/a/l 

l-g2 


(13) 


_ 1 +g^a^ 

'^XX ~ (l-g2)2 

C _ 

'^YY ~ (i-a2)2 


(14) 


Q_ — 

^XY ^YX 


(1-9^)' 


(“+ a) 


The minimal spectral density Smin corresponds to the quadrature = X cos(0min) + 

Y sin( 6 *miii), with 

1 C-... I c.... 

(16) 


^min = t : arctan 

Z (d) 


^XY + ^YX 


_ C-- 
W ^YY 


and is equal to 


5'min — 2 + ^YY 's/i^XX ^YY^ + i^XY + 

Using Eqs. (14), we obtain the explicit expressions for 6 *min and Smin- 


(16) 


^ TT 1 2 a 

c'min = w — V arctan ■ 


^(«2 _ 1) 


and 


(17) 




1 


7.2 


(1-^T 


1 + ^ + — 




'i: 


a^ 


a^ 


— +r\a+- 



(18) 


With our normalization, S'min is below its value with equal oscillators (i.e., S'min < 1/(1 + ^)^) 
whenever a 7 ^ 1, as shown in Fig. 4. We also underline in the same hgure that measuring 
the fluctuations in the antisymmetric quadrature (i.e., with 6 = x/4) is not the optimal 
choice whenever a 7 ^ 1. In our experiment, a < 1.1 and the difference between S'min and 
1/(1 + g)‘^ is small (e.g., 0.248 with respect to 0.25 for g 1), well below the experimental 
resolution. However, the accurate choice of 6 *min remains a crucial requirement even for such 
a weak asymmetry as soon as g >~ 0.75. 
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FIG. 4: Lower panel: with solid lines, we show the minimal spectral density S'min that can be 
obtained for the asymmetry parameter a = 1.1 (red), corresponding to our experimental param¬ 
eters, and with a = 2 (blue). The black dashed line shows S'min when a = 1 (i.e., for identical 
oscillators). With dashed-dotted lines, we report the fluctuations in the antisymmetric quadrature 
(with 6 = 7r/4). Upper panel: optimal phase 0min when a = 1.1 (red) and a = 2 (blue). 

Dynamic variance for the system of two different oscillators 


Below threshold, the spectra of the X and Y quadratures can be calculated from Eq. 8, 
translated into the Fourier space. Without parametric driving (with ^ = 0), the spectrum 
of X is 



Spi 1 1 

2 (miu;i)2^2+ 


(19) 
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and its variance is 


°° 0 ^ ^ da; 


K)= 




FI 


( 20 ) 


_ 2f 

For the Y qnadrature, the spectrum and the variance are given by Eqs. (19) and (20) by 
exchanging the subscripts 1 -H- 2. Before analyzing the effect of the parametric drive, we 
normalize the quadratures to the respective standard deviation calculated at ^ = 0. We thus 
dehne, similarly to what we have done for the spectral densities, X = Xj(j\ and Y = Y/ay- 
In the strongly quantum limit, cr^ = yjh/uimi and = ^h/u 2 m 2 : our normalization 
procedure corresponds to considering the usual dimensionless variables of the oscillators, 
normalized to the zero-point fluctuations. 

The evolution equation X + AX = ^ is now characterized by the system matrix 


A = 


71 

2 


^ga' 


f j/a' f 


( 21 ) 


where now the asymmetry parameter is 


a = 


SP2^l^l 


„ ( 22 ) 

OFim2UJ2 

that becomes a' ~ ^/Qi/Qi in the classical limit and a' ~ \/l 2 ph in the strongly quantum 
limit, and by the noise vector ^ = (y^/i, \^f 2 )- 

Below threshold, we can still calculate the stationary spectra 


= 

Fyp = 


+ I) 71 + I 


n;' + 4)l2 + #471 


/V 


/V 


S 


XY 


+ ^YX = -^4^^(7i«' + 72/a') 


(23) 


V = ( a ; 2 _ 7 m (1 _ ^ 2)^2 ^ 2 ( 71 ^ 


The maximally squeezed quadrature and the corresponding variance are derived by hrst 
calculating the spectral integrals cr~~ = (nnd similar expressions for a‘~~ and 

cr~~), then using Eqs. (15) and (16) where the spectral densities are replaced by the respec¬ 
tive covariance matrix elements. 

A different strategy for calculating the time evolving variance can be applied even above 
threshold, starting from the formal solution of the evolution equation, that reads 


X(t) = </>(t)X(0) + / (/)(t-r)^(r)dr 
Jo 


(24) 
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where 


e-^+* 0 


= M 


0 





(25) 


and M is the matrix of the eigenvectors of A, whose eigenvalues A± are given explicitly 
in Eq. (11). The elements of the covariance matrix for the vector X have the following 
temporal evolution 


= ^nW4x(0) + ^i2W4y(0) 



+7i / 0ii(^ - s)ds+ 72 / - s)<ls 


0iiW<^2i(^)4x(O) + <^12W<^22(^)4y(0) 


+7l / lAllC* - s)</'2i(< - s)ds+ 72 / I#>l2(i - s)i#>22(* - s)ds 


(26) 

(27) 


((t|;,-, and a~~ are found from the above equations exchanging X ^ Y and 1 eg 2), As 
in the stationary case, the maximally squeezed quadrature and the corresponding minimal 
variance are derived using Eqs. (15) and (16) where the spectral densities are replaced by 
the respective covariance matrix elements. 
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